The Universe contains approximately 6 times more dark 
matter than normal baryonic matter, and a directly observed 
fundamental difference between dark matter and baryons 
would both be significant for our understanding of dark mat- 
ter structures and provide us with information about the basic 
characteristics of the dark matter particle. We discuss one dis- 
tinctive feature of dark matter structures in equilibrium, namely 
the property that a local dark matter temperature may depend 
on direction. This is in stark contrast to baryonic gases. We 
used X-ray observations of two nearby, relaxed galaxy clusters, 
under the assumptions of hydrostatic equilibrium and identi- 
cal dark matter and gas temperatures in the outer cluster re- 
gion, to measure this dark matter temperature anisotropy /3dm, 
with non-parametric Monte Carlo methods. We find that /3dm is 
greater than the value predicted for baryonic gases, /3gas = 0, 
. at more than 3a confidence. The observed value of the tem- 
J> ' perature anisotropy is in fair agreement with the results of cos- 
O mological N-body simulations and shows that the equilibration 
of the dark matter particles is not governed by local point-like 
, interactions in contrast to baryonic gases. 
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Abstract 



1. Introduction 

The possibility that the temperature of the dark matter depends 
on direction is usually expressed through the velocity anisotro- 

py 
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where af and ct^ are the 1 -dimensional tange ntial and radial 
veloc ity dispersions in a spherical system (Binn ev & Tremaine 
19871) . Our intuition from classical gases leads us to loosely re- 
fer to these dispersions as the dark matter temperatures in the 
tangential and radial directions with respect to the centre of the 
galaxy cluster If most dark matter particles in an equilibrated 
structure were purely on radial orbits, then /3dm could be as 
large as 1 and, for mainly tangential orbits /3dm could be arbi- 
trarily large and negative. 

The temperature of a baryonic gas is only well-defined 
when the gas is locally in thermal equilibrium, and this gas 
equilibration is achieved through point-like interactions. In that 
case the gas temperature is independent of direction, which is 
expressed as /3gas = 0. This is because the mean free path, typi- 
cally in the tens of kpc dSarazinll 19861) . is much shorter than the 
cluster scale of Mpc. If dark matter was collisional and hence 
achieved equilibration through collisions, then one would also 
have /3dm = 0. 

Numerical N-body simulations of collision-less dark mat- 
ter particles show that the dark matter temperature anisotro- 
py is different from zero, and for galaxy clusters these sim- 
ulations show that /3dm goes from zer o in the central re- 
gion to 0.4-0.6 towards the outer regio n jCole&Lacevll 19961: 
Carlberg et al]ll997tlHansen & Moorell2006HHansen & Stadell 



20061 ). A mass-averaged /3dm is close to 0.3. 

During the assembly of galaxy clusters, the baryonic gas is 
shock-heated and eventually achieves energy equipartition with 
the gravitationally dominating dark matter at a temperature that 
is directly related to the temperature of the dark matter. This is 
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true as long as radiation and similar non-gravitational effects 
are negligible. The most sensible definition of dark matter tem- 
perature is by averaging over the 3 directions 
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and we describe in the next section how to derive this dark 
matter temperature from X-ray observations of the gas. Here 
ks is the Boltzmann constant, rup the proton mass, and fi the 
mean molecular weight (we assume /i = 0.61). We later dis- 
cuss to what extent the assumption of equality between gas 
and dark matter temperatures is supported by numerical sim- 
ulations (from an average over a large set of simulated clusters, 
some of which are significantly more perturbed than the ones 
considered here), and estimate the possible effect on /?dm- 

Our approach does not make any assumptions about the 
parametric form of ma ss or d ispersion profiles in contra st to 
earlier related studies ( Natar aian & Kneib 1996; Ikebe e Tal] 
|2004 . We later demand that the reconstructed dark matter tem- 
perature and the observed gas temperature must be equal in 
equilibrated regions that are not affected by radiative processes. 
This is a fair assumption since it only relies on the principle of 
equipartition between the dark matter and the gas. However, 
we emphasise that this is an assumption to be tested in the near 
future on high-resolution numerical simulations. 

2. Finding tlie darl< matter temperature 

For relaxed and spherically symmetric galaxy clusters, one can 
use X-ray observations of the hot, ionized gas to deduce the 
de-projected gas temperature and gas density as functions of 
radius. These are needed in the equation of hydrostatic equilib- 
rium (Fabricant et al. 1980; Sarazin 1986), which for spherical 
structures relates the total mass of gravitating matter at a given 
radius to the radial dependence of gas temperature and density 
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where G is the gravitational constant. We here assume that tur 
bulence is negligible, which will have to be te sted in the fu- 



ture e .g. by using line-broadening in metal lines dSunvaev et al 
2003h . 



We consider two highly relaxed clusters, which are likely 
to have only very little turbulence. That equilibrated struc- 
tures have reliably reconstructed mass (i.e. obey the hydro- 
static equilibrium) was support ed in the com parison between 
lensing and X-ray observations jAllen 1998 ). We also assume 
that non-gravitational entropy injection into the baryonic gas 
(e.g. from a central AGN) can be ignored in the region we con- 
sider Radio cavities are typically located within much smaller 
cluster-centric dist ances than those excised in our analysis (e.g. 



Birzan et al. 2004) We thus have both the total mass and the 



baryonic mass (since the cluster plas ma is optically t hin); since 
the stellar component is negligible ( Fukugita et al. 1 1998.) . we 
can easily get the dark matter mass and dark matter density as 
functions of radius. 

To make the connection with the d ark matter tempera- 
ture, w e must consider the Jeans equation dBinnev & Tremaine 
19871), which relates the dark matter density and velocity dis- 
persions with the total gravitating mass: 
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We are assuming that the Jeans equation for spherical struc- 
tures is accurate for dark matter, which has been shown to be a 
good assumption by numerical simulations dRasia et al ] |2004 . 
Equation (|4]i can be integrated to give 
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where p is defined from the dark matter density and anisotropy 
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Equations. (|3]|6ll provide all the tools needed to use the ob- 
served gas density and gas temperature to derive the dark mat- 
ter temperature. The only free parameter is the dark matter 
anisotropy, /3dm- We assume for simplicity that /?dm can be 
treated as a constant throughout the observed part of the galaxy 
cluster, a simplification that future data clearly will be able to 
lift. Thus, for any given value of /Sdm, we can calculate the DM 
temperature, which can then be compared with the gas tem- 
perature. This gives us the possibility of comparing different 
values of /?dm by standard statistical means. 

3. Two quiet clusters 

We used X-ray data from XMM-Newton of the two nearby 
galaxy clusters A 2052 and Sersic 159—03 (also known as 
Abell SllOl). These are highly relaxed clusters, where a tem- 
perature decrement in the central region is clearly identified 



in the de-projected data dKaastra et alj l2004t iPiffaretti et al 
20051) . These two clusters were chosen for three reasons. First 
of all, to extract the dark matter temperature non-parametiically 



one needs data at large radii. Secondly, the cluster surface- 
brightness map must be highly circular, with no evidence of a 
substructure, and the visual inspection of the de-projected tem- 
perature must show a smooth behaviour, indicating a fully re- 
laxed cluster. Finally, the de-projected gas temperature must be 
robustly observed. These clusters were analysed assuming th e 
standard cosmological ACDM model bv lPiffaretti et al.l d2005h . 
who show that the outer temperature decreases by approxi- 
mately 30% from the maximum temperature. These authors 
also show that a one-component gas temperature gives an ex- 
cellent fit to the observed spectra. Even though the gas temper- 
ature profiles of these two clusters are slightly different, they 
bot h appear to provide a fair fit to a universal temperature pro- 
file jPiffaretti et alJl2005h . 



4. Non-parametric analysis 

We treat the data in an entirely non-parametric way by Monte 
Carlo methods. The input data is a set of 7 radial values for 
the de-projected gas temperature and density and their corre- 
sponding error bars. These are determined from de-projected, 
spatially resolved spectra. A detailed description of the de- 
projection technique is given in lKaastra et al.l d2004l) . 

For each radial bin, we select randomly a temperature and 
density with a Gaussian-distributed value around the observed 
number and a width corresponding to the observed error bars. 
We then use the equation of hydrostatic equilibrium (Eq. [3]) to 
derive the total mass as function of radius. This total mass ex- 
tends only to a radius of 0.5-0.8 Mpc, so we select a random 
number between —4 and —2 for the logarithmic dark matter 
density slope at larger radii, which is a larger range than ex 



(5) pecta tions from both simulation and theory ("e.g. lDiemand et al 



20041) . The gas mass is already negligible beyond the outermost 
bin. We checked that varying these assumptions has virtually 
no influence on the results. 

We can now subtract the gas mass from the total mass. 
Having both the total mass and the dark matter density as func- 
tions of radius, as well as the assumed value for /?dm, allows 
us to integrate eq. Q and hence get the dark matter temper- 
ature from eq. (|2|i. At no point do we make any assumption 
about the form of the dark matter or gas profiles, nor about 
boun dary conditions. This is in contrast to earli er related stud- 
ies (.Nataraian & Kneiblll996t[lkebe et allllooi) . 

For each radial bin we now generate 10 000 models, and we 
can proceed with a frequentists statistical analysis. We calcu- 
late the median value of all the derived model temperatures, and 
we select the range covered by the central 70 — 75% of these 
models as representative of the error bars of the reconstructed 
dark matter temperature, ATdm- This number of central models 
is chosen to make the best per degree of freedom of order 
unity. We are being conservative since we are including more 
than the normal 68% of the models, corresponding to 1 stan- 
dard deviation for a gaussian probability distribution. As total 
error bar we use the quadratic sum of this reconstructed error 
and the observed temperature error bars, ATg (which is anyway 
much smaller). It is worthwhile to mention that a few of the 
produced models are non-physical in the sense that they may 
have decreasing mass locally; however, we choose the most 
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conservative approach and keep all models, hence slightly in- 
creasing the error bars while only very mildly shifting the me- 
dian. 

In Fig. 1 we show the observed gas and reconstructed dark 
matter temperatures for A 2052 (top-panel) and Sersic 159—03 
(bottom panel), with their corresponding error bars, for the case 
of Pdm = 0.6. We see from Fig. 1 that the gas and dark matter 
temperatures are in good agreement in the outer region. This 
value of (3 was chosen because it leads to reasonable agreement 
between the dark matter and gas temperatures. For a different 
value of /3, the reconstructed DM temperature will be different 
according to Eqs. (|2] |5]l. The main effect comes from Eq. |2] 
and shows that a smaller (or negative) (3 will imply a higher 
reconstructed temperature. 
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Figure 1 Observed baryonic (green stars) and reconstructed 
dark matter (red diamonds) temperatures with error bars for 
A 2052 (top) and Sersic 159—03 (bottom), as function of ra- 
dius. A dark matter velocity anisotropy of /3dm = 0.6 is as- 
sumed in both plots, because it gives a reasonable fit to the 
outer region. The radiative cooling is prominent in the central 
dense region, and in the analysis we include only the 4 outer- 
most temperature bins. 

There are large radial variations due to the non-parametric 
treatment; small non-monotonous valuations in the measured 



gas temperature and density imply non-trivial variation in the 
reconstructed dark matter temperature. Similarly, the error bars 
of the reconstructed dark matter temperature directly reflect the 
error bars of the measured gas density and especially tempera- 
ture. 

The central region of the cluster is very dense, and the gas 
temperature is mo st likel y governed by radiative and conduc- 
tive processes (e.g. lSarazin 1986) . so we expect the dark matter 
and gas temperatures to agree only in the outer part. The cor- 
responding radius where radiative processes become important 
is roughly where the temperature starts decreasing as one ap- 
proaches the centre. In particular, the central dark matter tem- 
peratures appear to have a quite different radial beha viour from 



the shape expected from numerical simulations (e.g. lEke et al 



1998c iRasia et al.ll2004h . This is clearly related to these being 
derived directly from the gas profile, which is strongly affected 
in the central region by non-gravitational processes. The most 
recent numerical simulations show good agreement between 
the simulated and observed gas temperatures in the outer re- 
gion ( Pratt et al. 2007 ). whereas the central region still allows 
for improvements in the simulations to reach agreement with 
observations. 

We have discussed how to reconstruct the dark matter tem- 
perature for a given assumed velocity anisotropy above. We 
therefore proceed and reconstruct the DM temperature for 
a range of different velocity anisotropies, and then compare 
them. 

We assume that gas and DM temperatures in the 4 outer- 
most radial bins must agree are do thus not include the region 
dominated by radiative processes. This allows us to perform a 

comparison between the reconstructed model and the ob- 
served temperature. We show the resulting figure of Ax^ as a 
function of /3dm in Fig. 2. One finds that a positive and non- 
zero dark matter anisotropy is preferred, with /3dm > 0.2 for 
Sersic 159-03 (and ^dm > for A 2052) at Ax^ = 9, coiTe- 
sponding to 3cr confidence. The two structures have /3dm < 1 
only at Ax^ = 3(7). This is the first time the dark matter tem- 
perature anisotropy has been measured, and it is comforting 
that it agrees fairly well with cosmological N-body simula- 
tions, which find /3dm ^ 0.3. We emphasise that we do not 
claim that the clusters have anisotropies of 0.45 or 0.7 (which 
naturally would disagree with numerical simulations), but are 
merely stating that the anisotropies are greater than zero and 
consistent with numerical simulations. 

For robustness we checked that keeping 68% of the models 
gives psr degree of freedom of 1.5 for A 2052 (and 1.1 
for Sersic 159—03), while keeping the best-fit points for /3dm 
virtually unchanged. In this case the resulting error bars of /3dm 
become slightly more stringent. 

We performed the same study using only the outermost 
3 bins in the statistical analysis, and the central values for 
/3dm were found to be near 0.4 for A 2052 and at 0.7 for 
Sersic 159—03, in excellent agreement with the results above. 
In this case one only keeps 40% of the models in order to get 
X^/d.o.f. close to unity, so we do not trust the resulting (slightly 
more restrictive) confidence interval for /3dm in that case. 
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Figure 2 Ax^ as a function of /3din- Both galaxy clusters have a 
preference for radial dark matter velocity anisotropy, with best- 
fit values at 0.45 (A 2052, red stars) and 0.7 (Sersic 159-03, 
blue diamonds). The 3a confidence level is shown with the hor- 
izontal (green dashed) line. A vanishing /3dm is ruled out at 
Ax^ of 9 and 18. 



5. Discussion 

The selected clusters appear to be fully relaxed; however, as 
noted above, the possibility of bulk motion in the gas still 
exists. Such a bulk motion implies that the reconstructed to- 
tal mass is underestimated since there will then be an extra 
term on the rh.s. of Eq. (3). Numerical simulations, includ- 
ing both gas and DM, have shown that for equilibrated struc- 
tures this underestimati on in the reconstructed mass is less than 
15% kasia et al.l2004 . This imphes that is underestimated 
by at most 15% (Eq. (5)). In addition, the dark matter temper- 
ature will be underestimated, since the energy of the gas also 
contains a term from the bulk motion. This underestimate is 



also at most 15% (iRasia et alJl2004l) in the region considered. 
These two effects partially (possibly almost totally) cancel 
out in the determination of the velocity anisotropy. Therefore 
/3dm is accurately determined, and the error from a possible 
bulk motion in the gas is significantly smaller than 20% (see 
Eq. (2)M 

We have here considered two clusters that both have a cen- 
tral temperature decrement, and we see that the reconstructed 
dark matter temperature is very different from the gas tempera- 
ture in the central region. There is growing evidence that clus- 
ters separate into two distinct classes, where one class has de- 
creasing central temperature, the so-called cool-core clusters 
like the two structures considered here, and the second class 
has a roughly constan t temperature in the central region (e.g. 
Sanderson et al. 20061) . It will be very interesting to perform 
another study similar to the one made here, to see if these non- 
cool-core clusters have dark matter temperatures that may even 
agree in the central region. 



' We are in the process of quantifying these estimates through high 
resolution numerical simulations, and the results will be presented 
elsewhere. 



It is worth mentioning that the temperature anisotropy can 
be measured in principle in an underground directional sensi- 
tive detector; howev er, it will require a lar ge dedicated experi- 
mental programme ( iHost & Hansenll2007b . 

Studies of hig hly non-equilibrated merging clusters (e.g. 
Clowe et alj|2006l) have shown that the dark matter has a much 
smaller scattering cross section than baryonic gas. We have 
here extended our understanding of dark matter to show that the 
equilibration of dark matter structures is not governed by point- 
like collisions. This demonstrates that dark matter behaviour is 
fundamentally different from baryons. This is particularly im- 
portant for understanding what drives structure formation and 
the evolution of cosmological structures. The standard theory 
of structure formation is based on the one basic assumption 
that dark matter is indeed collisionless, and we here provide 
observational evidence that this is a correct assumption. 
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